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Abstract 

Employing Maxwell's equations as the field theory of the photon, quantum 
mechanical operators for spin, chirality, helicity, velocity, momentum, energy and 
position are derived. The photon "Zitterbewegung" along helical paths is explored. 
The resulting non-commutative geometry of photon position and the quantum ver- 
sion of the Pythagorean theorem is discussed. The distance between two photons 
in a polarized beam of given helicity is shown to have a discrete spectrum. Such 
a spectrum should become manifest in measurements of two photon coincidence 
counts. The proposed experiment is briefly described. 



1 Introduction 

In recent years there has been considerable mathematical and physical interest in non- 
commutative geometry|l| and its physical consequences. Perhaps, the most simple 
example of a non-commutative geometry consists of the geometric plane described by 
Hermitian operator coordinates {X, Y) that do not commute; e.g. there exists an area 51 
such that 

i 51 (1) 

As a consequence of Eq.Q, there is a "quantization" of the Pythagorean theorem in 
that the distance D - VI 2 + P has a quantized spectrum|2|. For n — 0, 1,2,3, ... , 
we have 

(X 2 + Y 2 ) \n) = D 2 \n) where D„ = y/(2n + 1)51 . (2) 

Experimental examples of physical systems which can be usefully described by co- 
ordinates in the non-commutative plane include the following: (i) vortices in superfluid 
4 He films[ 3 J, (ii) quantum Hall effect charged magnetic vortices in two-dimensional 
electron liquids[4|, (iii) quantum interference phase between two alternative paths in 
the plane (as in the Aharonov-Bohm effect)|2| and (iv) high energy charged lepton 
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beams stored in cyclotrons||5]|6|[7). In all of these experimental examples, there ex- 
ists the ordinary Euclidean geometry of position vectors (i"i,r2, . . .) of the particles 
which make up the system. The Euclidean coordinates may appear in the quantum me- 
chanical wave functions of the system. Euclidean coordinates may also appear in the 
classical rulers employed to construct the lengths scales of the apparatus used in mea- 
suring the quantum mechanical system properties. However, the quantum coordinates 
(such as vortex positions or the positions of cyclotron orbit centers) are described by 
strictly non-commuting operators. Thus, both classical Euclidean geometry and quan- 
tum non-commutative geometry live tranquilly together in the theoretical description 
of the above experimental quantum systems. 

Such tranquility also exists in the quantum mechanical behavior of zero mass non- 
zero spin quantum mechanical particle beams. These include neutrinos, photons| 8 1 and 
gravitons. Our purpose is to discuss the case of zero mass photons ||9| ll0llTTlll2lll3l 
1141 . where the experimental technology 1 15 1 (of quantum optics) is most fully devel- 
oped. Our central theoretical results concern two photon coincidence counts measured 
from polarized beams. The statistics of the counts should reflect the non-commutative 
geometry of the photon positions. 

In Sec|2]we write Maxwell's equations in the vacuum in a form conducive to de- 
scribing the photon as a spin one object. The quantum mechanical spin, chirality, 
helicity, velocity, momentum and energy operators for a single photon will be derived. 
In Sec|3]the "velocity" operator of the photon will be explored. Writing the photon 
velocity operator as r = v = (v x ,v y ,v z ), one finds that each component formally has 
the possible eigenvalues -c, or + c. However, each component also has non-trivial 
commutation relations with the other components finally yielding 

v ■ v = v 2 + v 2 + v\ = 2c 2 > c 2 (photons). (3) 

Eq.0 indicates a considerable amount of "Jitterbugging" in the photon motion. For a 
Dirac particle, wherein Ypirac • ^Dirac = 3c 2 , Schrodinger called the resulting (superlu- 
minal) velocity fluctuations Zitterbewegung (German Jitterbugging). In more detail, it 
will be shown that the component of photon velocity perpendicular to the photon mo- 
mentum p undergoes a rotation at angular velocity Q. = c\p\/H. The resulting photon 
"path" may then be described as helical. Along the axis of the helix, the photon moves 
with velocity c. But the photon also moves around the helix in such a way that the total 
speed is greater than c. In Sec@J the Wigner photon position operator will be discussed. 
The photon position will be shown to obey the rules of a non-commutative geometry 
with an associated quantum Pythagorean theorem. The motivation for introducing such 
a position operator can be understood in terms of the photon angular momentum as is 
discussed in Sec[5] In Sec|6] the non-commutative geometry of a two photon system 
will be discussed. The distance between two photons in a polarized beam is shown to 
have a discrete spectrum. Such a spectrum should become manifest in measurements 
of two photon coincidence counts. In the concluding Sec0 the technology of these 
proposed experiments will be briefly described. 
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2 Maxwell's Equations 



The Maxwell equations describing a photon in the vacuum are well known. The rel- 
ativistic covariance is well established so that we pick a particular inertial frame (say 
the "laboratory frame") in which they read 

divE = 0, 

divB = 0, 

™ 1 / 5B \ 

curlE = — — , 

c \ at ) 



curlB 



c\dt 



Introducing the complex electromagnetic fields F + and F_ via 

f± = e + m, 

allows us to write Maxwell's equations in the more compact form 



div¥ ± = and i \ \ = +c curlF ± . 



(4) 



(5) 



(6) 



The operator curl defines the spin one matrices S = (S x ,S y ,S z ) according to the 
equivalent definitions 

U = curM 



and 
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Eq.(|8) implies the explicit representation 
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Note that the usual spin commutation relations, 

[S x , S y ] = iS z , [S y ,S z ] = iS x , [S Z ,S X ] = iS y , 

hold true as well as 

S 2 X +S* + S 2 Z =S(S + 1) where 5 = 1. 
Introducing the complex column three vectors F + and F_ via 



(E x ± iB x \ 

Ey + iBy 

E- + iB- 



(9) 



(10) 



(11) 



(12) 



3 



yields Maxwell's equations in the form 

i—± = +/c(S ■ grad)F ± (13) 
Ot 

where Eqs.(|6j, |7J and © have been invoked. Ea. (ll3> describes (in both righthanded 
F + and lefthanded F representations) the vacuum field equations for spin 5=1 mass- 
less photons. 

To see more clearly what is involved, let us consider the case of massless spin 
one-half particles, i.e. the Weyl massless neutrinos. Defining the momentum operator 

as 

p = -iti grad (14) 
and with cr = (<x v , <x v , cr.) denoting the Pauli matrices, 

'-(! i\ "-a o> h; .,)■ 

the massless "Weyl neutrino" equations read 

dil/+ 

in = c(cr ■ p)i// + righthanded Weyl, 

ot 

iti^- = -c{(T-vt)iIf- lefthanded Weyl. (16) 
ot 

The photon analogues to Eqs.(ll6> are 

dF 

iti—r— = c(S • p)F+ righthanded Maxwell, 
ot 

dF. 

in = -c(S-p)F_ lefthanded Maxwell. (17) 

Ot 

There is a difference between the spin one half Eqs. (I16> and the spin one Eas.( ll7> . 
A lefthanded and a righthanded neutrino are different particles. Only lefthanded neutri- 
nos have so far appeared in the laboratory which is a basis for parity violation. In spite 
of the wide spread theoretical righthanded cross product convention for writing down 
Maxwell's equations, electromagnetic theory is inherently parity symmetry invariant. 
Nevertheless, it is useful to construct a six component field column vector| 16 1 out of 
the three component righthanded and lefthanded column vectors defined in Ea.(ll2>: 

F = (£). (18) 

The operators for the photon are now in part described by 6 x 6 matrices naturally 
constructed by partitioned 3x3 sub-matrices. For example, employing the 3x3 spin 
matrices in Eq.(|5|l we may write for the 6x6 spin one matrices 



(19) 

The photon chirality matrix may be defined as 

HJ »). 
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The chiral spin £ matrix product 

« = Ts* = ( 8 _° s ), (2i) 

is useful for describing the massless photon Hamiltonian 

•K = c(a ■ p). (22) 

Eqs.jl7L dl 81 . fl21i and and ( 1221 imply the photon (Maxwell) wave equation in the 
Schrodinger-Dirac form 

dF 

ih— = <HF. (23) 
of 

To obtain the energy eigenvalues for the photon, one needs to solve 

W p (r) = e(p)F p (r). (24) 

Choosing /^(r) = F p {Q)e' p ' r ^ 1 yields the energy eigenvalues 

e+(p) = c|p| (photon), 
e_(p) = — c|p| (anti - photon), 

e (p) = (forbidden), (25) 

where the anti-photon (negative energy state) is merely the photon moving backward 
in time. For p + 0, the zero eigenvalue of energy is forbidden in virtue of the vacuum 
Gauss' law divE = and divB = 0. The photon velocity operator follows from the 
commutator 

▼ =I[W,r] = ^=ca = c(r 5 E). (26) 
« dp 

The photon helicity operator A is conventionally defined as 

A = fc^, (27, 



The allowed eigenvalues of helicity are A = ±1 while the eigenvalue A = is forbidden 
by the Gauss law constraints div¥ ± — of Eq.(|6}. The photon Hamiltonian is given by 

-H = (v • p) = c{a ■ p) = cr 5 (S ■ p), (28) 

which can then be expressed directly in terms of chirality Fg and helicity A as 

<H = c|p|r 5 A. (29) 



The eigenvalue spectrum of 'H in Eq.( I25> follows directly from the representation in 
Eq.@5)- 
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3 Photon Velocity 



Let us now contemplate how fast a photon is moving. For a photon of momentum p 
and energy e(p) = c|p| the mean photon (three vector) velocity is the expected light 
velocity value 

~WR de lp\ 
dp dp \\p\J 

On the other hand, from Ea. (l26> we find that photon velocity operator obeys 

v v = c 2 (r 5 i:) • (T S E) = c 2 (£ 2 + ly + E 2 ) = c 2 S(S + 1) (31) 

where [Fs, E] = and (Fs) 2 = 1 has been invoked. Since the photon has spin S = 1, it 
follows that 

v ■ v = 2c 2 . (32) 

Thus, the root mean square velocity of a photon is the superluminal value c y2. What 
is the nature of the velocity fluctuations? The answer resides in the operator for the 
photon acceleration a = v; i.e. 

v=-k"H,v] = kvp),v]. (33) 
n n 

Employing 

[v;,vy] = c 2 [r 5 S,-,r 5 S;] = c 2 [Z,-,E/] = ic 2 e ijk I, k = icT 5 e iJk v k (34) 

together with Eq.j33t yields the equation of motion for the photon velocity 

v = £2 x v wherein hSl = cF^p. (35) 

The physical meaning of Ea. (l35> is unambigous. The velocity of the photon precesses 
about the momentum p - hk direction at an angular velocity 

fi ± = ±(^) = ± c k (36) 

with the sign depending on the chirality Fs of the photon. 

From the classical electromagnetic theory viewpoint, one envisions either a right 
or left circularly polarized electromagnetic wave with frequency Q. = c|k|. From a 
photon viewpoint, with p = ftk, the photon moves with velocity c in the direction of k. 
The photon velocity component perpendicular to the k-direction is rotating (in either a 
right or left handed direction) at an angular velocity Q ± = +ck. The physical picture is 
depicted in Fig [2 

In position space, the photon moves on a "helical" path. The velocity component 
along the axis of the helix is light speed. On the other hand, the photon is going around 
and around on a helix with a total superluminal speed of c V2. Let us now consider in 
more detail how the helical path enters into the position representation for the photon. 
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velocity precession frequency Q=c|k| 

A 




Figure 1 : The circularly polarized photon velocity v precesses about the wave vector 
k = p/h axis at an angular velocity £1 ± = ±ck with the sign depending on chirality. 
The component of the photon velocity parallel to the wave vector k is light speed; i.e. 
c = (vk)/|k|. 



Photon Momentum p 
A 




Figure 2: For a photon of momentum p = ftk, we may point the momentum along the 
axis of a helix. The photon position r may than be written as a vector R pointing to the 
helix axis plus a vector u normal to the helix axis. The vector u moves in a circle at 
angular velocity il which depends upon the chirality via H£l = Tscp. 
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4 Photon Position and Helicity 



Consider the position of the photon as an operator in momentum space 

d 

r = ih—. (37) 
op 

The above Ea. P7> implies the usual Heisenberg commutation relations 

[Pk,rj\ = -iM Jk . (38) 

The position r of the photon can be decomposed into a sum of the position R of the axis 
of the helix and the displacement u from the helix axis to the helix coil. As depicted in 
FigE 

r = R + u, (39) 

where 

x p) _ (Mj\ vxp 
Ipl 2 ~ \ c ) |p| 2 

Radius vector u rotates in a circle remaining normal to the photon momentum p = fik 
according to the precession equation of motion 

li = -[7f,u] =ilxu wherein ft£l = cF 5 p. (41) 
n 

The radius squared of the helix is then given by 



» = ^=3^ = hr -=f • < 4C » 



For a spin 5 = 1 photon with helicity A = ±1, Eq. (l42> implies the square radius of the 
helix 

u u = |— j = A 2 (or equivalently) Vu u = |u| = A. (43) 

The radius of the helix orbit is completely determined by the wave length A = (2nX). 
However, the components of the radius vector from the helix axis to the helix path do 
not commute among themselves. In detail, Ea.J40> implies 

[Uj,ud = i(£-\ teju{r\\- (44) 



PI/ MP 

For example, if a photon moves with a fixed momentum p = (0, 0, ft /A) along the z-axis 
with helicity A = 1, then the components of u do not commute [u x , u y ] = iX 2 . Quantum 
mechanically, the radius of the helix is fixed at (u 2 . + u 2 ,) = A 2 but one cannot tell the 
values separately of the non-commuting components u = (u x , u y , 0). 

From Eqs.l|34}, OD. EDJ and ED it follows that [fl„ Uj\ = [n - u u Uj\ = 0, which 

together with [r,-, r^J = implies the non-commutative geometry of the coordinates R; 
i.e. 

"^-(afMS)- <45) 
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If a photon moves with a fixed momentum p = (0, 0, HI X) along the z-axis with helicity 
A = 1, then the components of R = (X, Y,Z) do not commute [X, Y] = -iX 2 and the 
quantum Pythagorean theorem of Eq.@ reads Vx 2 + Y 2 \n) - /tV2n + 1 \n) where 
n = 0,1,2,3 



5 Photon Angular Momentum 

The non-commutative geometry of the Wigner photon coordinate R can be motivated 
purely by considerations relating to the photon angular momentum J. As is usual one 
writes the total angular momentum as a sum of the orbital angular momentum and the 
spin angular momentum 

J = r x p + WL. (46) 

For a massless particle such as the photon, one seeks to express express the total angular 
momentum in terms of helicity. One employs the identity 

px(Exp) = |p| 2 {2I- Wa} (47) 



P, 

where the helicity A is defined in Ea.J27>. Eas. J46> and (1471 imply 

J = rxp-uxp + ^|^-)A. (48) 
\P/ 

where u is defined in Ea. (l40> . If we employ the Wigner photon position 

then the total angular momentum may be written in terms of the proper (zero mass 
particle) orbital angular momentum L = R x p and the helicity A according to 

J = L + fii:A = Rxp + /z(-^|A, 

IP/ 

A = +1 (allowed helicity), 

A = (forbidden helicity). (50) 
It is the Wigner photon coordinate which gives rise to a non-commutative geometry 

via 

[Rj,Rii] - -iMjk 

Mjk - (SfMwH^- (5i) 

The axial vector area, 

A=/ * 2 (^r) wherein lpl = ^' (52) 

plays an important role in the technology of photon detectors as will be discussed in 
what follows. The main physical point is that the position of a photon in the plain 
normal to the momentum is only defined within an area |A| - A 2 - |u| 2 . 
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photon 1 photon 2 




origin 



Figure 3: Two photons with the same momentum and helicity, move along different 
helical paths described by different Wigner coordinate vectors R\ and R2. The coordi- 
nate vector p — R2 - Ri describes the displacement vector from one helical axis to the 
other helical axis. 



6 Two Photon States 

Consider two photons each with the same momentum p and helicity A. Let Rj and R2 
represent, respectively, the Wigner positions of the first and second photon. Finally, let 



P = R 2 - Ri 



(53) 



represent the vector displacement from the helical axis of the first photon to the helical 
axis of the second photon as shown in Fig|3] From Eqs.( I51> . d52l > and J53i it follows 
that the displacement between helical axes obeys the non-commutative geometry rule 



[Pj,Pki = -iejki(An)i with area A i2 



' ft 2 PiA^ 
Pi 



r n 2 P2 A 2 ) 



pI 



(54) 



The distance between the two photon helical axes is then quantized in units of X ac- 
cording to 



Pi = P2 



= (0,0, ft/*), 

= 2A 2 (2n+ 1) |n> , 

= 0,1,2,3,... . 



(55) 



Now let us consider a monochromatic beam of photons of a given helicity inci- 
dent normal to a plane of photon detectors. Let us also consider two photon co- 
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incidence counts for the photon detectors in the the plane. According to our cen- 
tral non-commutative geometry Ea. J55l >. the Wigner distance between the two pho- 
ton coincidence counts should be quantized into spectrum D„ = A V2(2n + 1) where 
n = 0, 1,2, 3, . . . . Coincidence counting peaks at detector separations {D„} should be 
observable in the laboratory if the "pixel size" of the photon detectors are small on the 
scale of A = (A/2n) where A is the wavelength of the photons. 

7 Conclusion 

Rarely in experiments on optics does one try to resolve the photon position to within a 
spatial resolution smaller than A. For example, in the best commercial digital camera's, 
the pixel length scale L > A since experimentally it appears difficult to resolve photon 
positions on a smaller scale. To the authors' knowledge this lower bound on pixel 
resolution has not been attributed to the non-commutative geometry of photon positions 
but it does appear likely that such a geometry sets bounds on how well a photon can be 
localized. 

If the photon detectors are of the electric dipole type L « A, and if two photons in 
a circularly polarized beam are simultaneously measured, then the coincidence count 
distance D can in principle be measured to a resolution higher than A. Observation of 
peaks spaced at quantized distances D„ = A V2(2n + 1) would serve the useful purpose 
of corroborating the helix picture we have here theoretically obtained. 

The question then is: how easy would it be to obtain experimental evidence for 
the predicted effect? Common position-sensitive devices like charge conductive de- 
vices (CCD's) are not sensitive to single photons. Devices suitable for single photon 
detection based on photocathodes such as photomultipliers or hybrid photodiodes are 
unlikely to be fabricated with pixels small enough to lie inside a wavelength of light. 
Nevertheless, there is reason to be optimistic. 

Single photon detection requires a material with a bandgap E g less than he/ A. A 
useful number to keep in mind is that 1 eV corresponds to a wavelength of 1.2398 
microns, which is in the infrared. 

Fabrication requirements would suggest that it makes sense to go as far into the 
infrared as far as possible where there are essentially two options, both of which are still 
in rather early stages of development: 1) superconducting tunnel junctions (STJ's) ll7l 
or 2) pixellated avalanche photodiodes (APD's)| 18 1. 

Superconducting tunnel junctions have bandgaps of tens of meV, which correspond 
to breaking up Cooper pairs. They can, in principle, be sensitive to photons of a small 
fraction of an electron volt. This would allow pixels of several microns to be used, but 
to the best of our knowledge this has not yet been done. 

Avalanche photodiodes can be used for single photon counting and can also be fab- 
ricated as finely pixellated devices. Single photon counting has not yet been demon- 
strated with position sensitivity, but should certainly be feasible. Bandgaps are reason- 
able with common semiconductors like silicon and germanium having values less than 
an electron volt, and thus being sensitive to infrared of wavelengths longer than a mi- 
cron. Pixellation at the submicron level requires masks made using wavelength of light 
which are much shorter, but this is in principle possible with ultraviolet light, or with 
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shorter wavelengths becoming more and more available at synchrotron light sources. 

Quantum efficiencies for both STJ's and APD's can be quite high (in excess of 
60%), and noise should not be a problem if one uses them in a gated mode with a 
sufficiently intense correlated photon source. We hope to return to more detailed studies 
of the feasibility of an experiment along these lines in the near future. 
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